Stability of undercompressive shock profiles 



Peter HOWARD Kevin ZUMBRUN 

February 8, 2008 



Abstract 

Using a simplified pointwise iteration scheme, we establish nonlinear phase-asymptotic orbital stabil- 
ity of large-amplitude Lax, undercompressive, overcompressive, and mixed under-overcompressive type 
shock profiles of strictly parabolic systems of conservation laws with respect to initial perturbations 
|uo(a;)| < Eo(l + \x\)~ 3 ^ 2 in C 0+a , Eo sufficiently small, under the necessary conditions of spectral and 
hyperbolic stability together with transversality of the connecting profile. This completes the program 
initiated by Zumbrun and Howard in |53|. extending to the general undercompressive case results ob- 
tained for Lax and overcompressive shock profiles in |51|. |3U|. |53|. |55|. |44|. |37 | -|41 | . and for special 
undercompressive profiles in |31 | -|32 | . |23|. In particular, together with spectral results of :59 , our results 
yield nonlinear stability of large-amplitude undercompressive phase-transitional profiles near equilibrium 
of Slemrod's model |5L)| for van der Waal gas dynamics or elasticity with viscosity-capillarity. 

1 Introduction 

In the series of papers |55) . |37 j -|41 j . Zumbrun and Mascia-Zumbrun, building on methods introduced in 
EHl, have established nonlinear L 1 n H 3 — > LP (resp. L 1 n H 2 — ► L p ) orbital stability, p > 1, of large- 
amplitude Lax- type shock profiles of systems of conservation laws with viscosity (resp. relaxation), by a 
simple shock-tracking argument using mainly L q — > L p bounds on the linearized solution operator. For 
precursors of this method, see, e.g., JHj, [5H]-j57], |U _ E2j EH- EH- See also the alternative arguments 
carried out in [3U| for small-amplitude Lax- type shock profiles of systems with artificial (Laplacian) 
viscosity, and in JS3j , 03] for large-amplitude Lax- or overcompressive- type profiles of systems with general, 
possibly degenerate viscosity. 

The purpose of the present work is to point out that a simple pointwise version of the argument of [55] , 
|37 j -|41 | may be applied also to under-, over-, and mixed under-overcompressive shock profiles of strictly 
parabolic systems, giving a simple and unified treatment of shock stability independent of the amplitude or 
type of the connecting profile, depending only on the necessary, Evans- function condition established in |53j . 
|14| . 0], equivalent to spectral and hyperbolic stability plus transversality of the connecting profile 

as a solution of the associated traveling- wave ordinary differential equation (ODE). In particular, we obtain 
for the first time nonlinear stability of general undercompressive profiles such as arise in phase-transitional 
gas dynamics and elasticity 0B]-[IS] or multiphase flow E3> extending results obtained for special 

undercompressive profiles in [23] 

Moreover, the slight additional detail afforded by our pointwise description is sufficient to give also 
convergence of the phase shift, or phase-asymptotic orbital stability (definition recalled below), which was 
lacking in [S3], [37|-02]. This completes the one-dimensional stability program initiated in [53], at least for 
systems with strictly parabolic viscosity, giving a complete characterization of stability analogous to that 
obtained by Sattinger in the scalar case |45| . The assumption of strict parabolicity is physically appropriate 
for applications to phase-transitional shock waves, at least in one dimension; see Remark^] Indeed, together 
with the spectral analysis of [3JJ, our results yield one-dimensional nonlinear stability of large-amplitude 
phase-transitional profiles near equilibrium of Slemrod's model for van der Waals gas dynamics or elasticity 
with viscosity-capillarity, one of only two large-amplitude stability results that have so far been obtained 
for physical models (the other being stability of profiles of isentropic gamma-law gas dynamics with 7 = 1 
|3ti| ) . Stability of undercompressive profiles for systems with degenerate viscosity remains an interesting 
open problem. 
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Consider a traveling-wave solution 



u(x, t) 



u{x — st), 



z — >±oo 



lim u(z) = u± 



(1) 



or "shock profile" , of a system of conservation laws 



ut + f(u) 



x 



(B(u)u x ) 



(2) 



x, teM.,u,fe K", b e w 



corresponding to an "ideal" , or discontinuous shock wave 



u(x, t) 



u_ x < st, 
u + x > st, 



(3) 



of the associated hyperbolic system 



Ut + f(u)x = 0. 



(4) 



Without loss of generality (changing to coordinates moving with the shock), take s — 0, so that (JJJ becomes 
a stationary, or standing-wave solution convenient for stability analysis. 
Following |S7], we make the standard assumptions: 

(HO) f,BeC 2 . 

(HI) Rea{B)>0. 

(H2) &(f'(u±)) real, distinct, and nonzero. 

(H3) Re a{-ikf(u±) - k 2 B{u±)) < -9k 2 for all real k, some 6 > 0. 

(H4) There exists a solution u of C3>~©, nearby which the set of all solutions connecting the same 
values u± forms a smooth manifold {u s }, S <G U C u° — u. 

Definition 1. An ideal shock @ is classified as undercompressive, Lax, or over compressive type according 
as i — n is less than, equal to, or greater than 1, where i, denoting the sum of the dimensions £_ and i+ of 
the center-unstable subspace of rf/(u_) and the center-stable subspace of df(u+), represents the total number 
of characteristics incoming to the shock. 

A viscous profile is classified as pure undercompressive type if the associated ideal shock is under- 
compressive and I = \, pure Lax type if the corresponding ideal shock is Lax type and £ = i — n, and pure 
overcompressive type if if the corresponding ideal shock is overcompressive and £ = i — n, I as in (H4) and i 
as in Definition^ Otherwise it is classified as mixed under-overcompressive type; see t<5ffi . 

Pure Lax type profiles are the most common type arising in standard gas dynamics, while pure over- 
and undercompressive type profiles arise in magnetohydrodynamics (MHD) and phase-transitional models. 
Mixed under-overcompressive profiles are also possible, as described in [^21, [SSj, but seldom encountered; 
indeed, we do not know a physical example. In the pure Lax or undercompressive case, {u s } = {u(- — 5)} 
is just the set of all translates of the base profile u, whereas in other cases it involves also deformations of 
u. For further discussion of existence, structure, and classification of viscous profiles, see, e.g., [32], |53j . 



as defined in (H4), for any solution u of 0) with initial data sufficiently close in some norm to the original 
profile u. If, also, the phase 5(t) converges to a limiting value 5(+oo), the profile is said to be nonlinearly 
phase-asymptotically orbitally stable. 

An important result of |53| was the identification of the following stability criterion equivalent to L 1 —> L' p 
linearized orbital stability of the profile, p > 1, where D(X) as described in denotes the Evans 

function associated with the linearized operator L about the profile: an analytic function analogous to the 
characteristic polynomial of a finite-dimensional operator, whose zeroes away from the essential spectrum 
agree in location and multiplicity with the eigenvalues of L. 
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(2?) There exist precisely I zeroes of D(-) in the nonstable half-plane Re A > 0, necessarily at the origin 
A = 0. 

As discussed, e.g., in under assumptions (H0)-(H4), (V) is equivalent to (i) strong spectral 

stability, <r(L) c {Re A < 0} U {0}, (ii) hyperbolic stability of the associated ideal shock, and (iii) transver- 
sality of u as a solution of the connection problem in the associated traveling-wave ODE, where hyperbolic 
stability is defined for Lax and undercompressive shocks by the Lopatinski condition of |33|-|35|. jl()j and 
for overcompressive shocks by an analogous long- wave stability condition 56]. Here and elsewhere a denotes 
spectrum of a linearized operator or matrix. 

The stability condition holds always for small-amplitude Lax profiles J7|, GH|> EB1> EH]- |^4| . |43|. 
but may fail for large-amplitude, or nonclassical over- or undercompressive profiles pQ, ^2], E0|> EH- ^ 
may be readily checked numerically, as described, e.g., in 0, El- It was shown by various techniques 

in [27|-0T], 031 that the linearized stability condition (T>) is also sufficient for nonlinear orbital 

stability of Lax or overcompressive profiles of arbitrary amplitude. However, up to now, this result had not 
been verified in the undercompressive case. 

In this paper, we present a simple pointwise argument applicable to shocks of any type, establishing that 
(T>) is sufficient for nonlinear phase-asymptotic orbital stability. More precisely, denoting by 

af < 4 < • • • < at (5) 

the eigenvalues of the limiting convection matrices A± :— df(u±), define 

6{x, t) := C 1 + t)-V2 e -l-°7*l 2 /^ + ( J + t)- 1 /2 e -l-4*l 2 /" i ( 6 ) 

a~<0 a+>0 



M*, t) •■= X(x, t)Y(l+ \x\ + t)- l ' 2 (l + \x- aTt\y 1 ' 2 



A^i v 1 1 ^ \ 1 3 

Q-<0 

X(x,t) J2 (l + \x\+t)-^(l + \ X -a+t\ 

a+>0 



-1/2 



(7) 



(8) 



and 

i> 2 {x, t) := (1 - X (x, t))(l + \x ~ <hA + t^ 2 y 3/2 

+ (l-x(x,t))(l + \x-a+t\+t 1 / 2 )- 3 / 2 , 

where x(x,t) — 1 for x £ [a^t, a+t] and zero otherwise, and L > is a sufficiently large constant. 
Then, we have the following main theorem. 

Theorem 1. Assuming (HO)-(H^) and the linear stability condition (D), the profile u is nonlinearly phase- 
asymptotically orbitally stable with respect to C 0+a initial perturbations \uq(x)\ < Eq{\ + |a;|)~ 3 / 2 , Eq suffi- 
ciently small. More precisely, there exist S(-) and S(+oo) such that 

\u(x, t) - u s ^ {x)\ < CE {8 +^ + i/t2)(x, t), 

\S(t)\<CE (l+t)-\ (9) 

|<y(t) - <y(+c»)| < CEvii + ty 1 ' 2 , 

where u denotes the solution of (0 with initial data uq = u + uq. 

In particular, Thcorcm^yields the desired result of nonlinear stability in the undercompressive or mixed 
case, effectively completing the one-dimensional stability analysis initiated in |53|. 

Remark 1. Pointwise bound @ yields as a corollary the sharp L p decay rate 

\u{x,t) -u s{t \x)\ LP <CE (l + tyi {1 -p\ l<p<oo. (10) 
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Remark 2. The profile 9 may be recognized as the superposition of Gaussian "approximate diffusion waves" 
moving along outgoing characteristic directions, while the profiles ip\ and "02 respectively account for nonlinear 
interactions occurring within the characteristic cone [a^t, a^t] and the algebraically decaying tail of the initial 
data. The profiles 9 and ipj correspond roughly to the nonlinear diffusion and linearly coupled waves 9 and rj 1 
in the more detailed description of the solution carried out for initial data with the same decay rate in ]30f . 
\5'i$ for the Lax and overcompressive case. The latter works estimate u — u, which contains an additional 
error 

u(x) - u m (x) ~ (du s /dS)S(t) = 0(e-^){l + ty 1/2 ) (11) 

near the shock layer that is not present in our analysis. 

A difference of the undercompressive from the Lax or overcompressive cases is that the time-asymptotic 
distribution of mass is no longer determined in a simple way by the mass of the initial data, making difficult 
the description of nonlinear diffusion and coupled waves. We avoid this difficulty by estimating only joint 
upper bounds and not the size or shape of component waves. Besides undercompressive stability, this yields 
also considerable simplification in the pointwise analysis of Lax and overcompressive profiles. In particular, 
we nowhere attempt to identify cancellation in our estimates of nonlinear interactions, taking into account 
only transversality of interacting signals, similarly as in (different-family) Glimm interaction estimates for 
the hyperbolic case \1 ftf . Compare with the analyses of 1511 1.Vfll \44\ in which a crucial aspect is to identify 
cancellation in the computation of the linearly coupled wave r\. Compensating for the lower resolution in 
our scheme at the level of diffusion waves is the higher resolution afforded by shock tracking, as reflected in 
the absence of term <|11[1 . That is, by sufficiently resolving the nondecaying lowest-order part of the Green 
function corresponding to shift in the shock location, we are able to ignore the details of higher- order parts. 

Remark 3. Multidimensional nonlinear L 1 nL°° — > L p stability, p>2, has been established using L q — > L p 
resolvent bounds for Lax and overcompressive shocks in all dimensions d > 2 but for undercompressive 
shocks only in dimensions d > 4; see 156$ - 15S$ . As discussed in 121$ . }5fif . \22j . stability of undercompressive 
shock fronts in physical dimensions d — 2 and 3 remains an open question even in the (diffusive-dispersive or 
diffusive-higher order diffusive) scalar case, for which detailed pointwise Green function bounds are available. 
An interesting future direction might be to attack this problem by pointwise methods similar to those of this 
paper. 

Remark 4. In one dimension, a change of coordinates reduces Slemrod's model for van der Waals gas 
dynamics or elasticity with viscosity- capillarity to a 2 x 2 parabolic system with diagonal, strictly parabolic 
viscosity, to which our results may be applied; see ]5(A \5iJ$ . Likewise, in one dimension, the standard models 
for imiscible three-phase flow in porous media may be expressed as a 2 x 2 parabolic system with viscosity 
that is strictly parabolic on the interior of the physical state space (the set of saturations summing to one) 
and degenerate on the boundary; see '31. Thus, our results generically apply here, too. However, note that 
shocks with one state on the boundary do arise in Riemann problems of physical interest and are not covered 
by our theory, nor is the degeneracy of the symmetric constant-multiplicity type encountered in real viscosity 
models. This would be an interesting case for further investigation. In multi- dimensions, neither of these 
transformations is possible, and so a more special analysis of each specific equation would be required for a 
stability analysis. 

Plan of the paper. In Section we recall the linearized estimates carried out in Assuming 
certain pointwise convolution estimates, we carry out in Section|31the nonlinear stability analysis of the Lax 
and undercompressive case. By a slight modification of the argument, we carry out in SectionQ]the nonlinear 
stability analysis of the complementary Lax and overcompressive case. In Sectional we establish stability of 
mixed-type shocks in the special case B = constant , and in Scction|S|in the general case. Finally, in Section 
we carry out the deferred convolution estimates, completing the analysis. 



1 Defined in 1511 but not explicitly mentioned in either reference, r] refers to the contribution of quadratic source terms 
involving 9 alone, i.e., the second iterate in the nonlinear iteration through Duhamel's formula. 
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2 Linearized estimates 

We begin by recalling the pointwise linearized estimates established in [53| . |39j . expressed in a streamlined 
form convenient for the nonlinear analysis to follow. Linearizing J2J about a fixed stationary solution u s *(-) 
gives 

u t = L s *u:={B s *u x ) x -{A 5 *u) x , (12) 

where 

A s *u:=df(u s *(x))u-dB(u s *(x))(u,u s x *), B 5 * :— B{u 5 * (x)). (13) 

Proposition 1 ([53j, [39j). Under the assumptions of Theorem^ the Green function G{x,t;y) associated 
with the linearized equations (I12|) may be decomposed as G = E + G, where 

£ (^) = E^^ ejfot), (14) 
j=l i \s=s, 

ej (y,t) = J2 [errfn I V -±M ) - errfn ( V -^M ) ) l Jk (y) (15) 

for y < and symmetrically for y > 0, with 

\1%\<C, \(d/dy)lf k \<C 7 e-^, (16) 

and 

\dl y G{x,t-y)\ < 

n 

C(t~^/ 2 + K| 7 e-*l + \a x \e-^){^^ 1/2 ^ (x ^ a " t)2/Mt ^ VX+ 

fe=i 

i V V f -\/2 -{x-aJ(t-\y/a-\)f/Mt - V x+ (17) 



a, >0,a. <0 



, V v .-1/2 -(x-a+(t-|a/a^|)) 2 /Mt -jjx"^ 

^ X {\%t\>\y\} u e J' 

a^>0,af>0 

< I CK j < 2 for y < and symmetrically for y > 0, /or some 77, C, M > 0, where are as m Theorem^ 
> 0, x^ denotes the positive/negative part of x, indicator function Xs\ a -t\>\y\] * s ^ -^ or l a fc^l — If I an< ^ ^ 
otherwise, and 7 = 1 in t/ie mixed or undercompressive case and m i/ie pure Lax or overcompressive case. 
Moreover, all estimates are uniform in the supressed parameter S r . 

Remark 5. We will refer to the three differently scaled diffusion kernels in jJ7| ) respectively as the convection 
kernel, the reflection kernel, and the transmission kernel. The x-derivative estimate will only be required in 
the case of mixed-type profiles (see Section 6). Finally, we recall the notation 

1 f z 2 
errfn(z) := — / e~^ df . 

W ^ Loo 



Proof of Proposition ^ This is a restatement of the bounds established in 1211 for pure un- 
dercompressive, Lax, or overcompressive type profiles; the same argument applies also in the mixed under- 
overcompressive case. Also, though it was not explicitly stated, uniformity with respect to 5* is a straight- 
forward consequence of the argument . | 
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Remark 6. From (|f 5fl and <|R)I) . we obtain by straightforward calculation (see Jffffi) the bounds 
M».t)l<C£ Urr/n [ ) - errfn / 1 " °*' 



a, >0 



40* */ V\/W 



l e j(y 5 *) - ej(y,+oo)| < Cerrfn( ^— somea > 



e -|*+o* t r/ Ar ' 



o t >0 

e~ 



a, >0 



+ Cje l,y ' errfn — — err/n 1 



(18) 



/4 W vv 4/3 ^ 

a->0 

|^(y,t)|<C(i- 1 +7t- 1 / 2 e-"^l) £ e -\y+a^f/Mt 

a->0 

for y < 0, and symmetrically for y > 0, where 7 as above is one /or undcrcompressive profiles and zero 
otherwise. 

Remark 7. The main difference between the estimates for the mixed or undcrcompressive case 7 = 1 and 
the pure Lax or overcompressive case 7 = is the presence of slower-decaying e _e ' y ' terms in derivative 
estimates for e, G. As discussed in f^j ^ -fA 1 ^ . f5ffi. ■ ]55f . these are not only technical artifacts, but 
reflect real differences in behavior in the undcrcompressive case: specifically, that shock dynamics are not 
governed solely by conservation of mass, as in the Lax or overcompressive case, but by more complicated 
dynamics of front interaction as indicated by rapidly decaying modes ~ e~ e, ' y '. 

3 Stability of Lax or undercompressive profiles 

We now carry out the proof of Theorem ^ in the Lax or undercompressive case, which may be treated by a 
particularly simple argument. In these cases £ = 1, and u s — u(x — S), so that we may conveniently work 
with the "centered" perturbation variable 

u(x,t) :=u(x + S(t),t) -u(x), (19) 

for which becomes 

u t - Lu = Q(u, u x ) x + 8(t)(u x + u x ), (20) 

L := L°, where 

Q(u,u x ) = 0{\u\ 2 + \u\\u x \) 

Q(u,u x ) x = 0(\u\\u x \ + \u x \ 2 + \u\\u xx \) 

so long as \u\ remains bounded. 

Recalling the standard fact that v! is a stationary solution of the linearized equations Q12[l. Lu' — 0, or 



, : , , (21) 



G(x,t;y)u x (y)dy = e u x (x) = u\x), 



we have by Duhamel's principle: 

/>OC 

u(x,i) = G(x,t;y)u (y)dy 



v 

J-00 



G y {x,t— s; y)(Q(u,u x ) + Su)(y, s) dyds + 8{t)u'(x) 
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Defining 

S(t) = - / e{y,t)u Q {y)dy 



ft /- + oo 

(22) 

+ 11 e y (y,t - s)(Q(u,u x ) + 5u)(y,s)dyds, 

JO J-oo 

following [S3], [37|-[3H|, where e is defined as in ((T3|) (that is, e = J2j e j)i an d recalling the decomposition 
G = E + G, we obtain finally the reduced equations: 

/oc 
G(x, t;y)u Q (y) dy 

rt 'IS, (23) 



G y (x, t - s; y){Q[u, u x ) + Su)(y, s)dyds, 

and, differentiating (1221) with respect to t, and observing that e y (y,s) —r as s — > 0, as the difference of 
approaching heat kernels: 



£(*) = - / e t (y,t)u (y)dy 

*t ~?+oo _ (24) 

+ e yt (y,t- s)(Q(u,u x ) + 5u)(y,s)dyds. 



./-oc 



We shall make use of the following three technical lemmas, the proofs of which are given in Section 

Lemma 1 (Short-time theory). Under the assumptions of Theorem^ for data uo G C 0+a (x), equations 
(12211 H23J) (alternatively, H5Q(I . (|4f \ of the following section) admit a unique local solution u G C 0+a (x) n 
C 0+a / 2 (f) 7 5 G C* 1+a / 2 (i), extending so long as \u\ c o+ a remains bounded. Moreover, on this domain, 
sup 2 + "01 + ^2) (■£,•) remains continuous so long as it and \S(1 + t)\ are uniformly bounded and, 
for t > t > sufficiently small, 

sup \u x \(9 + i>! + ^ 2 )^ 1 (z,<) < Ct- 1 / 2 sup \u\{9 + Vi + ^ 2 y\z, t - r). (25) 

Z Z 

Lemma 2 (Linear estimates). Under the assumptions of Theorem^ 

f +00 

\G(x,t; y)\(l + |y|)- 3 / 2 dy < C(0 + Vi + V> 2 )(x, t), 

|e t (j/,t)l(i + |y|)- 3/2 dy<C(i + t)- 3 / 2 , 

«+oo 

\e{y,t)\{l + \y\)-" 2 dy<C, 
\e(y,t) - e(y, +oo)|(l + \y\V 3/2 dy < C(l + ty 1 / 2 , 
for < t < +00 , any a, M > 0, for some C > 0, where G and e are defined as in Proposition^] 



(26) 



7 



P. HOWARD and K. ZUMBRUN, Stability of Undcrcompressive Shock Profiles 



Lemma 3 (Nonlinear estimates). Under the assumptions of Theorem^ 

ft /- + OC 

G v (x,t- s;y)\V(y,s) dyds < C{6 + fa + ip 2 )(x, t), 

^-00 

t />+oo 

/ \e yt (y, t - s)\*(y, s) dyds < C(l + t)~\ 

J -co 
+oo />+oo 



r-t-oo r-tcc 

/ / \e y (y,+o )\V(y,s)dyds<C 7 , (27) 



t />+oo 

/ \e y (y, t-s)- e y (y, +oo)|tf (y, s) dyds < C(l + £)- 1/2 , 

J-oo 



oo p 4- oo 



/or 



K(y,£- s)|¥(y,«)dy <C(f +t)- 1 / 2 , 



*(y, s) := (1 + s) l / 2 s- l / 2 {6 + ^ + ^) 2 (y, *) 
+ (1 + s)- 1 (0 + ^ 1 +7A 2 )(y, S ) 



(28) 



Proof of Theorem ^ Lax or undercompressive case. Define 

C(t):= sup (lup + ^+^r^O + I^Kl + s)). (29) 

y,0<s<t K ' 

We shall establish: 

Claim. For all £ > for which a solution exists with £ uniformly bounded by some fixed, sufficiently 
small constant, there holds 

m <c 2 (E +at) 2 ). (30) 

From this result, provided Eo < 1/46*2, we have that £(£) < 2C 2 Eq implies £(£) < 2C2E0, and so we may 
conclude by continuous induction that 

C(t) < 2C 2 E Q (31) 

for all £ > 0. (By Lemma ^ u G C 1 exists and £ remains continuous so long as ( remains bounded by some 
uniform constant, hence (|3f I) is an open condition.) Thus, it remains only to establish the claim above. 

Proof of Claim. We must show that u(9 + ipi + V^) -1 and | j(s)|(l + s) are each bounded by C(Eq + ^(£) 2 ), 
for some C > 0, all < s < t, so long as C remains sufficiently small. 
By combined with we have for t > 1 that 

\u x (x, £)| < Cat - + Vi + ih)(x, t - 1) 

<C 2 C(£)('5 + ^i+V'2)(a;,£) 1 j 

and for < t < 1 that 

Ms, t)l < Ct-VWW + Vi + lka)(a?, 0) 
<C 2 C(£)£- 1/2 ('5 + Vi+V'2)(x,£). 
Combining these estimates, and recalling definition (|29|l . we obtain for all £ > and some C > that 

i^)i<cw(i+r\ 

KM)| <C(*)(fl + 4>i+^2)(x, £), (34) 
K(z, £)| < CC(£)(1 + £) 1/2 £- 1/2 (0 + ^1 + fo)(x, £) 

and therefore 

|(Q(«,u x )+j«)(tf,«)|<C*(tf,«) (35) 
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with 4" as defined in (JS8J. 

Combining (|35fl with representations (|23[) ~H24 [) and applying Lemmas [21 and |3 we obtain 



\u(x,t)\ < I \G(x,t;y)\\u (y)\dy 

ft poo 



J-oo 



\G y (x,t- s;y)\\(Q(u,u x ) + 5u)(y, s)\dy ds 
<E I \G(x,t;y)\(l + \y\y 3 / 2 dy 



+ CC(<) 2 / / \G y (x,t-s;y)\*(y,s)dyds 

JO J-oo 

<C(So+CW 2 )(^ + ^l+V'2)(x,t) 



and, similarly, 



>Vi / \e t (y,t)\\u (y)\dy 

ft /* + OG 



\e y t(y,t- s)\\(Q(u,u x ) + Su)(y,s)\ dyds 

'0 J-oo 

pt p + oo 

<l \e t (y,t)\(l + \y\)- 3 / 2 dy+ / \e yt (y, t - s)\V(y, s) dy ds 



J-oo 



<c(^o + CW 2 )(i + _1 - 

Dividing by (8 + ipi + ip2){x,t) and (1 + respectively, we obtain l|30[) as claimed. 

From 1)30(1 . we obtain global existence, with ((t) < 2CE$. From the latter bound and the definition of ( 
in (|29Jl we obtain the first two bounds of 0. ft remains to establish the third bound, expressing convergence 
of phase 5 to a limiting value <5(+oo). 

By Lemmas El together with the previously obtained bounds (|35[) and C < CEo, and the definition (|29|) 
of C, the formal limit 



oc 

oo />+oc 



8(+oo) := / e(y,+oo)u Q (y)dy 

e y (y, +oo)(Q(u, u x ) + Su)(y, s) dyds 
; / E \e(y,+^)\(l + \y\)-V 2 dy 

C-Eo|e y (2/, +oo)|*(y, s) dyds 



OC' p-+-oo 



< CE 

is well-defined, as the sum of absolutely convergent integrals. 
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Applying Lemmas [2HB1 a final time, we obtain 

/oo 
|e(y,t)-e(»,+oo)||uo(i/)|di/ 
-oo 

ft r + oo 

I \ e y(y,t- s ) - ey(y,+oo)\\(Q(u,u x ) + Su)(y, s)\ dy ds 

J -oo 



oo /> + oo 



\e y (y, +oo)||(Q(u, u x ) + 5u)(y, s)\ dyds 
: / \e(y,t)-e(y,+™)\(l + \y\y 3 / 2 dy 

\e y (y,t- s) - e y (y, +oc)\$ (y,s)dyds 

J-oo 

oo r+oo 



\e y (y,+oo)\ty(y,s)dyds 

<c£ (i+tr 1/2 , 

establishing the remaining bound and completing the proof. a 



4 Overcompressive profiles. 

We may treat the overcompressive case by a slight modification of the argument of Section [3J which applies 
also to the Lax case. As the Lax and overcompressive case have already been treated by different means in 
|53j . [H| we shall only sketch the changes necessary for the argument, omitting most details. 

Modified equations. In the overcompressive case, I > 1, u s consists not only of translates of u, but 
also of orbits distinct from u. In particular, the different representatives are not all derived from a group 
action, and so we cannot use a centering transformation as in p9|l . which consists of the group operations 

Tg(u - T_$u) = T s u - u, 

where T a v(x, t) := v(x + ct,t) denotes translation in x. Accordingly, we work with the primitive variable 

u(x,t) :=u(x,t) -u m (x) (36) 

and center the equations instead, about some strategically chosen <5* , obtaining in place of lj20|l the modified 
perturbation equation 

u t -L 5 'u = Q s * (u, u x ) x + 8(t){du s /dS)\ S , + R 5 ' (S, u, u x ) x + S 5 ' (8, S t ), (37) 

where Q * is as in J2H) and 

R 5 ' = (A{u s *(x)) -A{u s ^{x)))u+ (b(u s '{x)) ~B(u m {x)))u x 

^ ' ^ ' (38) 

= ( e -n\y\\6-5*\(\u\ + \u x \)) 

and 

S s * =S((du 5 /d6) m -(du 5 /dS) lSt ) =0(e-^\5\\8~5*\) (39) 
account for "centering errors"; see, e.g., [T^], |2T], [S7| for related computations. 
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Defining 



rt r + oo 

6(t):=6,- I e(y,t)u 5 '(y)dy- / e(y, t - s)S 5 *(S, S t )(y, s)dyds 

-oo JO J — oo 

t p+oo 

e y (y, t - s)(Q s * (u, u x ) + i? 5 * (6, u, u x ))(y, s)dyds, 

J -oo 
S, ~ -8, 

u := uq — u 



(40) 



OO 

t poo 



by analogy with we obtain reduced equations 

POO 

u(x,t) = / G(x,t;y)u s Q *(y)dy 

G(x 1 t~ S ;y)S 5 *(6,6 t ))(y, S )dyds (41) 
G y (x,t- s;y)(Q s *(u,u x ) + R s * (S,u,u x ))(y, s)dy ds 

JO J -oo 

and 



./-oo 
t poo 



6(t) = - / e t (y,t)u 5 '(y)dy 



t r +oo 

e t (y,t- S )S 5 *(5,S t )(y, S )dyd S (42) 



J -oo 
t r + oo 

/ e yt (y,t- s)(Q s *(u,u x ) + R s * (8,u,u x ))(y, s) dy ds. 

J -oo 

Asymptotic shock location. An important feature of the Lax or overcompressive case is that the 
stability criterion (T>) in these cases implies that the (nonlinear) L 1 - asymptotic state of the perturbed shock 
must be formally determined by conservation of mass, in the sense that relation 

(u (y) - u(y))dy = (u +°° (y) - u{y))dy 

a+>0 aj<0 

is full rank, hence uniquely soluble for , S +00 for f^(u — u)(y)dy sufficiently small, where and 

denote eigenvalues and right eigenvectors of A± — f'(u±), denotes asymptotic mass in the jth outgoing 
characteristic field at ±oo, and <5 +oc denotes the asymptotic shock location, or, equivalently, 



/+oo 
(«o " u s ^)(y) dy= m t r t ( 44 ) 



for S +OQ = O(E ), the latter estimate a consequence of full rank; for further discussion, see |33 _ IS2) 1531, 
|56j . and references therein. Centering about 5* = 5+oo, we may thus arrange that 



/+oo r+oo 
u 5 '(y) = / (u - u 5 *)(y) dy=Y, m t r t> ( 45 ) 



while maintaining our assumptions on initial perturbation u * . In these coordinates, we may expect that 
\S(t) — <5»| decays to zero. 
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A second consequence of (1431) . this time at the linearized level, is that e(y, +00) (constant in the Lax or 
overcompressive case; see Proposition ^| must be orthogonal to all "outgoing modes" r+, > and and 
rj, a J < 0, hence with choice of coordinates (|45J) we have 

/•4- 00 

e(y,+oo)u 5 -(y)dy = 0. (46) 
A final consequence of 143|) is that the map 

-00 

U(u s - u)(y)dy € M £ (47) 



must be invertible, where II € I x n is any (constant) full rank matrix with rows orthogonal to outgoing 
modes . Reparametrizing by (5 = 5, we may thus arrange that 8 = f_™ TL(u s — u)(y) dy, and thus 



j 

U(du 5 /5)(y)dy = I e (48) 
for S in a neighborhood of the origin, so that f_ °° US S * (y : s) dy = for all s, and therefore 

r +00 

efo,+oo)S*'(y,»)dtf = 0. (49) 
Combining (|46fl and (|49|) . we obtain the alternative representation 

/oo 
(e(y, t) - e(y, +oo))u (y) dy 
-OO 

- / / (e(y,t- s) - e(y,+oo))S s *(S,8 t )(y,s)dyds (50) 

JO J-co 
pt p + oc 

+ 11 e y (y,t- s)(Q 5 *(u,u x )+R s *(8,u,u x ))(y,s)dyds, 

JO J-oo 

from which we may observe decay in \5 — 5*| without a priori knowledge of the global behavior of u. 

Stability argument. Working within the framework of equations (|41JI . (|42[) . and H50I) . we can carry out 
the proof of Theorem 0f or the Lax or overcompressive case by essentially the same argument presented for 
the Lax and undercompressive case in the previous section using Lemmas H |H1 together with the following 
Lemma proved in Section 

Lemma 4 (Auxiliary estimates). Under the assumptions of Theorem^ 

ft r + oo 

\G y (x,t - s;y)|$i(y, s) dyds < 0(9 + ^! + i/j 2 ){x,t), 

t /> + OC 



y 

J-oo 



J-oo 
t r +oc 



\e yt (y,t-s)\$x(y,s)dyds < C(l + f)~\ (51) 
\e y (y, t - s)|$i(y, a) dyds < C(l + ty 1 ' 2 



y 

J-oo 

t r + oo 

/ \G(x, t - s; y)\$ 2 (y, «) dyds < C(6 + ^ + ^ 2 )(x, i), 

J-oo 

£ /> + oo 

/ |e i (y,t- S )|<i> 2 ( 2/ , S )dyd S <C(l + t)- 3 / 2 , (52) 

J-oo 
t r+oo 

/ |e(y, t - s) - e(y, +oo)|$ 2 (l/, *) dyda < C(l + t)~ 3/2 , 

J-oo 
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where 

$ x (y, a) := 8-^(0 + + tfa)(y, *) < Ce-^s^il + s)-\ 
^ 2 (y, S )~e-^(l + s)-^ 2 . ( ' , ' >1 



(54) 



Specifically, defining 

C(i) := sup + ^ + ^ 2 )- 1 (y,t) + |j(s)|(l + s) + \5(s) - *„|(1 + s) 1/2 ), 

V,0<s<t K ' 

noting that 

|Q tf «|<C 2 *, < C 2 (* + *i), |S 5 -|<C 2 $2, 

and applying our convolution lemmas, we may obtain l|3U|) as before, yielding at once global existence and 
the claimed rates of decay. 

Remark 8. It was crucial in the argument to linearize about the limiting profile u s +°° in order that error S 
be manageable, i.e., 6 — > <5* as t — ► +oo. 

5 Mixed type profiles, constant viscosity case. 

We now present an alternative proof subsuming Lax, undercompressive, overcompressive, and even mixed 
under-overcompressive cases in a single argument. For clarity of exposition, we first restrict to the simpler 
case B = constant , which permits also the following slightly stronger result. The general case is treated in 
Section 

Theorem 2. Let B = constant . Then, assuming (H0)-(H4), and stability condition (D), the profile u is 
nonlinearly phase-asymptotically orbitally stable with respect to (not necessarly Holder continuous) initial 
perturbations \uq(x)\ < E Q (1 + |:e|)~ 3 / 2 , Eq sufficiently small. More precisely, © is satisfied for some S(-), 
<5(+oo), where it denotes the solution of with initial data Uq = u + u . 

Proof. Defining u as in i|36|) . we obtain 

u t - L s -u = Q s * (u) x + S(t)(du s /dS) ls , + R s ' (5, u) x + S S ' (5, S t ), (55) 

where 

Q 5 '=0(\u\ 2 ), (56) 

R s * = (a(u s *(x))- A{u 5 ^ (x))) u = 0{e-^\ \5 - 6*\\u\) (57) 
and S s * is as in (I39[l . Defining 

r t r +oo 



8(t)=6*- I e{y,t)u *(y)dy - / e(y,t - s)S * (5, S t ){y, s)dyds 

J — OO JQ J — OG 

+ / / e y (y,t-s)(Q 5 *(u)+R s '(5,u))(y,s)dyds, 

JO J -oo 



(58) 



u s *(y) := (u-u s *)(y), 



we obtain 



u(x,t) = T u (u,5,6,5*)(t) := I G(x, t; y)u * (y) dy 

J — OO 

/ G(x,t- S ;y)S 5 *(6 7 S t ))(y,s)dyd S (59) 

J — oo 

/ G v (x, t - s; y)(Q 5 ' (u) + R s * (S, u))(y, s)dy ds 

J-oo 
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and 

/oo 
£t{y,t)u s * {y)dy 
-oo 

rt p + oo 

/ et(y,t- s)S 5 '(S,St)(y,s)dyds 

JO J -oo 



e yt (y, t - s){Q 6 ' (u) + R 6 ' (S, u)){y, s) dy ds. 

^-00 



Defining now 

T 5 (u,6,5,5*)(t) := 



and 



(e(y, t) - e(y, +oo))u S ' (y) dy 

pt p + oo 

/ / (e(y,i- s) - e(y, +oo))S 5 "(S, S)(y, s)dyds 

JO J -oo 
rt p + oo 

/ / (e y (y,t-s)-e y (y,+^))(Q 5 '(u) + R 6 *(S,u)){y,s)dyds 

JO J -oo 

-oo />+oo 

/ e(y,+oo)S s -(S,5)(y,s)dyds, 

J — OO 
+ 00 /* + OG 

e y (y, +oo)((Q 5 * (u) + R K {5, u))(y, s) dy ds, 
^5+oo('M' , M*)(*) := 

f°° -8 

8*- e(y,+OQ))u *(y)dy 

J — OO 

00 p-\-oo 

/ e(y, +oo)S 5 "(5, S)(y, s)dyds 
J -00 

-00 r+oc 

/ e y (y,+oo)(Q 5 '(u) + R 5 * (S,u))(y, s) dy ds, 
J — 00 

we may express the solution of (158(1 - 1)60(1 equivalently as the solution of the fixed-point equation 

(u,S,S) = {T u ,T s ,Tg)(u,6,5,5*). 

in combination with 
Defining norm 

\(f,g, h)\ c ■= 1/(0 + i'i + M^U-c-.t) + \g(t)(i + *) 1/2 U» ( t) 
+ IM*)(i + *)U-(t) 

and Banach space 

B:={(/,fl,/»): !/,<?,% <+co}, 
we find by the estimates of the previous sections that, for 

\u -u\(x) <E Q (l + \x\)-^ 2 , 

Eq sufficiently small, (T u ,Ts,Tg,Ts +oo ) is a well-defined mapping from 

B(0,r) cBxI^BxM 

for r > sufficiently small, with 

\T\ BxM = O(E + \5*\ + \(u,6,6)\l). 
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Moreover, essentially the same estimates yield that T is Frechet differentiable on 5(0, r), with 



88 

and 



= 0(\(u,6,S\ c ), (68) 
o(u, 0, 0) 

0(1), (69) 



(71) 



-J±2. = 0(\(u,6,8\ c ), (70) 
where the final equality JJJJ follows from 

/oo 
e(y, +oo))(as 5 /^)| 5t (y) dy + 0(\(u, 8, 6\ c ) 
-oo 

= h + 0(\(u,S,5\ c ). 

In turn, relation 

/>oo 

e(y,+oo)){du s /d5)\ 5 ._{y) dy = h 

— OO 

follows from the standard fact that L 5 * (du 5 / d5)\$_^ — 0, hence 

+00 

G(x, t; y)(du s /dd) lSt (y) dy = (du 5 /dd) lSt (x) 

-00 

which, together with the fact that E = {du 5 / d8)\$^ (x)e(y, t)) represents the only nondecaying part of 
G(x,t;y) under stability criterion (T>), yields 

/+00 
e(x,+oc;y)(du 5 /d5) { M dy = (du 5 /d5) {s ,(x) 
-00 

in the limit as t — > +00. 

Combining (|68|I ~H70 |) . we find that, for Eo, r sufficiently small, T is contractive with respect to norm 

\(i,j,k,l)U :=\{i,j,k)\ c + C\l\, (72) 

for C > sufficiently large, with |T(0)|„ = O(E ). Applying the Contraction Mapping Theorem, we find 
that have a unique solution in B, from which the stated decay estimates follow by definition of 

Mm 



Remark 9. The above contraction mapping argument may be recognized as an alternative version of the 
Implicit Function Theorem argument commonly used to establish orbital stability, as for example in \2U$ . 

M 



6 Mixed type profiles, general case. 

The constant-viscosity assumption of the previous section made it possible to work in a weighted L°° norm, 
since we needed to gain only a single derivative in the associated nonlinear iteration scheme. To treat 
the general case, we work instead in a Holder space using Schauder-type smoothing estimates like those of 
Lemma ^ 
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Proof of Theorem ^ general case. In the general, variable-viscosity case, we may express the 
perturbation u of l|36|) again as the solution of equation |JSSJ, but with Q s * and R s * now depending also on 
u x . 

To accomodate this fact, we impose on uo also Holder continuity, 

K| c o+« < C, (73) 

and impose on u the uniform bounds 

\u\ c °+« < C, \u\ c i < Cr? +a , \u\ C 2 < Ct- 1+a , (74) 
or equivalently |tt| Q < C, for 

|u| a := sup (\u\ c0+a + \u\ c is?- a + \u\c2S 1 - 01 ) , (75) 

s>0 v ' 

together with \(u,S,5)\^ i: < E sufficiently small, for 

+ \djti (1 + t)-i{6 + + ^r'lL-i^t) (76) 
+ \g(t)(l + i) 1/2 Uo. (t) + \h(t)(l + t)\ L oo W 
(note: now augmented with derivative bound). Denote 

B x := {(f,9,h) : \f,g,h\ Cl < +00} 
C:={(f,g,h): \f\ a < +00}. 



(77) 



We now define Tg, 7], and 7s ^ by (|60() - (|62|l . exactly as before. However, we define T u now implicitly, as 
the solution of 

/oo 
G(x, t; y)u s * (y)dy 
-00 

- / / G(x,t- S ;y)S s *(6,6 t ))(y, S )dyd S 
Jo J-00 (78) 

ft pOO 

/ G y {x,t-s;y)(Q s '(u,d x T u ) 
J-00 

+ R 5 ' (8, u, d x %))(y, s)dy ds, 

or equivalently as the solution v of 

vt - L s ' (u)v = Q s * (u) x + 5(t)(du s /dS) lSt + R s * (6, u) x 
+ S S * (S,S t ),+f s *(u,u x ,v x ), 



(79) 



(81) 



where 

L s * (u)v := B{u s * + u)v xx - (A(u s * )v) x (80) 

and 

f 5 '(u,u x ,v x ) = 0((\u\ + \u x \)\v x \), 
T s *(u,u x ,v x ) x = 0((\u x \ + \u xx \)\v x \ + (\u x \ + \u xx \)\v xx \) 
for \u\ sufficiently small. 

Observing that L(u) may be expanded as a nondivergence-form operator with C 0+a coefficients, we may 
obtain by standard Schauder or parametrix theory, similarly as in Lemma ^ both short-time existence for 
(|79|l and also the smoothing estimates 

\v\c(t) < Ct- 1+a (\u \ c ° + « + \u\ a ), 

x (\u\ Cl + sup \u \(8 +Vi +V^) _1 (^,0)) 
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for < t < 1 and 

|«|c»(*) < C(\u \ c o +a + \u\ a + \v\ c i(t- 1)), 

sup\v x \(9 + i; 1 +yj 2 r 1 (z,t)<C 
z [S6) 

x (ju| Cl +supH(0 + Vi +^ 2 ) _1 (2,i- 1)) 

for f > 1. Substituting these bounds in (|78|l . and performing estimates similarly as in Section we find by 
the same type of continuous induction scheme that T is well-defined, and bounded from a sufficiently large 
ball about the origin to itself in {B\ flC) x IR . Moreover, similar estimates yield that T restricted to this 
same ball in (B\ n C) x IR 1 intersected with a sufficiently small ball in Boo x IR is contractive in the rescaled 
Qi norm 

\(i,j,k,l)U*:=\(i,j,k)\ Cl +C\l\, ( 84 ) 

with |'7"(0)|b oo xir < CEq. Combining these facts, we obtain a unique fixed-point solution in (B\ H C) x IR 1 , 
for which 

\(u,S,S)\ c <CE , 

yielding similarly as in Section [S] a global solution of the perturbation equations satisfying the claimed 
bounds. We omit the details. ■ 



7 Technical Lemmas 

We complete our analysis, finally, by the proof of the deferred technical lemmas used in Sections 
7.1 Short-time theory 

Proof of Lemma ^ We carry out the proof for equations i|22[l - 1[23[) . The proof for equations JSUJ, 
(|41Jl goes similarly. By standard Schauder or parametrix theory [151 1 X r> j . there exists a unique solution 
u e C 2+a (x) n C 1+ ^ (t) for t > sufficiently small of the original (unshifted) equation J2J), extending so long 
as |u|c°+°(a;) remains bounded and satisfiying uniform bounds 

\u\ C o + a (x) < C, \u\ cl{x) <c( T ^) _5+Q , \u\ C 2 {x) <c( T ^) _1+Q (85) 

depending only sup lul^o+o.^). 

This in turn determines 5, hence 5, through l[24[l by a straightforward contraction-mapping/continuation 
argument. (Note: in establishing contractivity for small time of the righthand side of (|24[l . we must establish 
bounds of form 

rt r + co _ i-t 

C \e y \{\u x \ 2 + \u\\u X x\) dy ds < C (\u x \ 2 L „ + \u\ L °° \u xx \ L ™) ds < 1 

JO J-oo JO 

for t > sufficiently small, which follow by \e y \i,i < C, a consequence of Remark [S] and integrability of 
l^zl^oo + (wie^Iloo , in turn a consequence of estimates From l)24[l . moreover, we easily obtain that 

8 € C 1+ ~ , and uniformly bounded for short time, so that perturbation u{x + S(t), t) — u(x), or equivalently 
its shift u(x,t) — u(x — S(t)), enjoys the same regularity properties as u. (Recall that u € C 3 , as a solution 
of the C 2 traveling-wave ODE.) Moreover, |u| c o+c, remains uniformly bounded so long as |u|co+ Q does, and 
vice versa, since their difference u(x — S(t)) is uniformly bounded in C 0+Q . 

This verifies the first assertion. To verify the second, observe for fixed to , r that u(x, t) for to — r < t < to 
by Duhamel's principle satisfies 

/oo 
G(x, t - (t - r);y)u(y, t - r) dy 
-OO 

+ / G y (x,t - s;y)(Q(u,u x ) + Su)(y,s)dyds 

JtQ—T J — OO 
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and therefore 

u x (x,i) = I G x (x,t - (t - T);y)u(y,t - t ) dy 



G x (x,t- s;y)(Q(u,u x ) x + 5u x )(y, s) dy ds. 

to — r J —oo 

Combining (|21(1 and (|74|) . we obtain 

pOO 

\u x (x,t)\< / \G x (x,t- (t - T);y)\\u(y,t -r)\dy 



+ C / \G x (x,t-s;y)\ ( 86 ) 

J t(j — T J — OO 

x (\u x \(s - (t a - r))-i+ a + \u\(s - (t Q - r))- 1+Q ) (y, s) dyds, 
from which we readily obtain 

\u x V- x (x,t)\ <(t- (t a - r))"5 aup\u(z,t - t)*- 1 ^, i - t)| 

z 

for to — r 5: t < to by a contraction argument based on the convolution estimates of Lemmas EH3 thus 
verifying (|25|) . (By uniqueness, the fixed point obtained by contraction mapping must in fact be u.) See 
Lemma 5.1 | 52| or Lemma 11.5, |53| for similar arguments. 

Likewise, we obtain by quite similar argument the Holder bound 

\df u\\^-\x, t)\<(t- (to - r))- a sup \u(z,t - t)*- 1 ^, to - r)|, 

z 

where d t 2 u(x, t) denotes limsup e _^ \u(x, -)| c o+f j t _ e t+ , with some abuse of notation. Combining this fact 

with the uniform bound l 1 ^ 1 ]/ -1 ! < C obtainable by direct calculation, we obtain the claimed continuity of 
sup 2 \u^ _1 (z, -)|. Indeed, we obtain Holder continuity, C 0+ ^. a 



7.2 Integral estimates 

Throughout the analysis, we will make use of the following lemmas. 

Lemma 5. Let f(y) > be a bounded, nonincreasing function on R + , and also let f G L^R). Then for 
any a > and z > 0, and for any u> > 1, 

X3 r— 

a^e-W f(y)dy < (\f(z/u)) A (a l ' 2 \\f\\ L , (K) ) 
o A 

+ [(^ll/llL~( R ))A(a 1 /2|| / || il(R) )]e-^ 2 , 

for any 7 < (1 — ^) 2 , and where A represents minimum. 

Proof of Lemma [SJ Lemma [5] is proven as Lemma 6.3 in t 

The following two lemmas, useful in analyzing the 9(x,t) terms in the nonlinearity 'J, can be proven in 
straightforward fashion by completing an appropriate square. 

Lemma 6. For any x, y, s, t, Mi, M 2 a, and b, we have 

(x — y - a(t - s)) 2 (y — bs) 2 (x — a(t — s) — bs) 2 
Mi(t - s) ' M 2 s ~ Mi(t-s) + M 2 s 
Mi(t - s) + M 2 s ( (xM 2 s - (aM 2 + bAh)(t - s)s) \ 2 



/ (xm 2 s - (aM 2 + oivii){t - s)s)y 
M 1 M 2 s(t-s) V M x (t - s)+M 2 s / 
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Lemma 7. For any x, y, s, t, M\, M 2 a, b, and c we have 

(x-%y-a(t-s)) 2 , (y~cs) 2 _ {x - a{t - s) - c±sf 



M t (t - s) M 2 s Mi(* - s) + M 2 (f ) 2 s 

. 6 ' 



Af 1 (i- S ) + M 2 (f) 2 s 



M 1 M 2 (f) 2 s (t-,s) 
(o(f ) 2 M 2 + (o(t - s) + cf s)M 1 )s(t -a)- x(f) 2 M 2 s^ 2 



1 b V Ml (t-s) + (%) 2 M 2 s 

Proof of Lemma El In each case of Lemma EJ we proceed in the case x, y < 0. The case y < < x is 
similar to the reflection estimate for x,y < 0. The case y > is entirely symmetric. (See Remark [3] for a 
discussion of our terminology regarding convection, reflection, and transmission kernels. We will designate 
the case 7 = as the Lax case and the case 7 = 1 as the undercompressive case.) 

For the first estimate in Lemma we consider the integrals 



\G(x,t;y)\(l + \y\y 3 / 2 dy. 
Convection estimate. For the convection kernel, according to Lemma we have 

t-We- 1 "*"*"**' (1 + \y\r 3/2 dy 

< C{t-^ 2 A (1 + \x - al t|)- 3 / 2 + (1 + t)- 1/2 e" i= #) , 

where M' > M can be taken as close to M as we choose (by choosing the 7 of Lemma sufficiently close 
to 1). In the event that x and have opposite signs, we have decay of the form (|a;| + t)~ 3 / 2 , which 
can be absorbed by the claimed estimates. In the event that x and have the same sign, the terms 

(1 + t)~ 1 / 2 e «i are exactly the 0(x, t). In order to see that the expression t~ x / 2 A (1 + \x — a^t\)~ 3 / 2 
can be absorbed into the sum 9 + tpi + ip 2 , we first observe that for \x — aZt\ < Cy/i, 

r 1 ' 2 A (1 + \x - a-t\)- 3 ' 2 < d(l + ty^e- 1 -^ 1 , 
for some constant C\. On the other hand, for \x — aZt\ > Cy/i, 

(l + \x- a- <|)- 3 / 2 < t-^ 2 (l + \x- a-t\)-^ 2 , 

which is sufficient in the case \x\ < |a7|t. Finally, for x S 2a^t], i -1 / 2 < C(|a;| + t)~ 1 / 2 , while for 
x (jL [-^-t, 2a^t], there can be only limited cancellation between x and a^t, and we have 

(1 + \x - a-t\)- 3 / 2 < C(l + \x\ + ty 3 ' 2 . 

For |x| > \&i\t, we analyze the case \x — a^t\ < Cy/i as above, leaving the case x < a^t — Cy/i, for which 

r 1 ' 2 A (1 + \x - a^t\y 3 ' 2 < C(l + \x- aZt\ + t 1 ' 2 )- 3 ' 2 . 

Reflection estimate. In the case of the reflection kernel, Lemma [S] provides the estimate 



< 



;-i/2 e %j (1 + \ y \y 3 / 2 dy 

c^r 1 ' 2 a (1 + \x - ay\y 3 ' 2 + (i + ty^e-'^^) , 
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which can be analyzed exactly as above. For the transmission kernel (and x < 0), we have exponential decay 
in x. In the event that \x\ > et for some e > 0, we have exponential decay in both t and x, which can be 
subsumed. In the event that \x\ < et, we have 

f -i/2 e k — t e -nl*l(i + \y\)- 3 / 2 dy 

OO 

which can be subsumed into the claimed estimates. The cases for x > are similar. 
For the second estimate of Lemma |5J we consider integrals of the form 

f \e t (y,t)\(l + \y\)- 3 / 2 dy, 

J — OO 

where according to Remark [5] 

\e t {y,t)\<Ct- 1 / 2 Yl e 

Q->0 

According to Lemma [5] we can estimate this integral as, 

\e t (y,t)\(l + \y\)-^ 2 dy < C [ r 1 /^-^^ (i + \y\y^dy 

J — OO 

< c(i + t)- 3/2 . 

The third estimate of Lemma |21 follows directly from the boundedness of e(y,t). 
For the final estimate of Lemma [3J we consider integrals of the form 

f | e ( y ,t)- e ( y ,+oo)|(l + |y|)- 3 / 2 , 

J —OO 

where according to Remark El 

\y\ - at 



\e(y, t) — e(y, +oo)| < Cerrfn 
for some C, M, a > 0. Computing directly, we have, then 



My/t 



|e(y,t)-e( 2 /,+oo)|(l + |y|)- 3 / 2 dy 



/ + OC 
errfn 
-OO 



i r° , r=£3r 
c 

2tt 



oo «/ — OO 



I f 2 / f MVt 

= —c 



2tt 



OO 





e- z2 dz)(l + \y\)- 3 / 2 dy 

t 

( / ^ e-* 2 dz)(l + \y\)-V 2 dy 



-y-at 



This completes the proof for x, y < 0. As discussed above, the cases in which either x or y is positive 
follow similarly. | 
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Proof of Lemma |3J For Lemma [31 the proof of each estimate requires the analysis of several cases. 
We proceed by carrying out detailed calculations in the most delicate cases and sufficing to indicate the 
appropriate arguments in the others. In particular, we will always consider the case x,y < 0. The case 
y < < x is similar (though certainly not identical) to the reflection case for x, y < 0. The estimates for 
y > are entirely symmetric. For the nonlincarity W we observe the inequality 

*<C(l+ S ) 1 / 2 .S- 1 /2 (0 2 + ^2 +v ,2 ) ^ 

Nonlincarity 9 2 . We begin by estimating convolutions of the form 

,-t /-O 

\G y (x, t-s;y)\(l + s) 1 ' 2 ' S - l 'H{y, sfdyds. 



IQ J-oc 

Lax convection, x < 0. For a~ < 0, we consider integrals of the form 

t /*0 (x — y — a~(t — s)) 2 (y — a~ s) 2 

/ (t - s) _1 e mij^s) (! + s )-V2 s -i/2 e xk—dyds, 

J-oo 

where the constant M arising in 9 2 is larger than the exact constant, L/2 (making the term an upper 
estimate). According to Lemma [SJ we have 



J -co 



(x- H -a fc (t-s))^ (y-a. s) z 

(t - s)~ x e (1 + s)- 1/2 s- 1/2 e T^dyds 

(t - a)- 1 e- l ""''^T' 7 ' i \ l + S )-V2 S -V2 



t r 



(x-a k (t-a) 

(t-s)- L e- 

J-oo 

t . x a -(q-+a-) S (t-a) 

x e ~ Ms(t-,) (v —t ) dyds 



< cr 1 ' 2 [\t - sy^ii + g )-v2 e - ( *~ a * (t i/ t ) " a ' 3>2 ds . 

Jo 



The convection rate can be either positive or negative, so we divide the analysis into three cases: (1) 
aj < < , (2) < aj < 0, and (3) aj < < 0. For the first, we observe that for \x\ > \aj\t, we can 
write 

x-a^(t- s) - ajs = (x- ajt) - (a^ - aj)(t - s). 

In this last expression, x — ajt and — {aj — aj)(t — s) are both negative and cannot cancel, and we can 
compute 

cr 1 ' 2 [\t - s)- 1/2 (i + s)" 1 /^- — 'm?'" 3 ' 3 ' 2 ^ 
t 







KCr 1 ' 2 / ( t - S )-V2( 1 + s) -l/2 e 







1/2 



< C{l + t)- l ' 2 e 

(We observe that the seeming blow-up of t~ x / 2 as t — ► is compensated for by the interval of integration 
s G [0,i].) In the event that \x\ < \aj\t, we further divide the analysis into cases s £ [0,t/2] and s £ [t/2,t). 
For s € [0, i/2], we write 

x — a^(t — s) — ajs = (x - ajt) + (aj - aj)s, 

for which x — ajt < and (aj — o,J)s have opposite signs and cancellation occurs. In this way, we have the 
balance estimate 

(x-a-(t- S )-a- S ) 2 

(l + s)-^ 2 e m 3 

f , ^-a7{t-3)-aJ s) 2 (a.--a-t) 2 

< C (1 + \x - ajt\)-^ 2 e + (1 + s)-V3 e tpt- 
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We observe that this kind of balance estimate is contained within the proof of LemmaEl but that the lemma 
cannot quite be directly applied, because we first must recognize which variables balance. In this way, the 
current analysis is a refinement of the analysis of We also observe that it is critical that for s £ [0, t/2] 
the cancellation comes from s (which on this interval does not yield t decay), while for s £ [t/2,t] it will be 
critical that the cancellation comes from (t — s). Indeed, at a purely technical level, this observation drives 
the analysis. Computing directly, we now have 

rt/2 (x-a7 (t-s)-a~ s) 2 

t- 1 ' 2 (t- s )-V2 (1 + s) -i/2 e fc m > ds 



< cr 1 ' 2 f'\t - s)- 1 / 2 [(l + \ X - dj;t\)- 1 '*e- l '~° hl *™~' 1 ' )a da 



(l + s)- 1/2 e 



ds. 



Finally, observing that in this case |x — a k t\ > c(\x\ + 1), and that t > |x|/|a ■ | we determine an estimate by 

which can be subsumed into V'l- I n the case s £ [t/2, t] we argue similarly, beginning with the relation 

x - a k (t ~ s ) - a J s = (x- ajt) - (Ofc - aj)(t - s). 

For the second case (aj < a^ < 0), we again begin with the subcase |x| > \aj\t, for which the argument 
for the case aj < < holds. Similarly, for the subcase \x\ < \dj~\t we write 

x ~ a k (* - s ) - a J s = (x- alt) + (aj^ - aj)s, 

for which again there is no cancellation, and we obtain an estimate by 9. In the event that \a~l\t < |x| < \a~ \t, 
we proceed in the cases s £ [Q, t/2] and s £ [t/2, i] exactly as in the case a J < < . The third case, 

< aj < follows similarly, though here we begin with the cases |x| > \aj~\t and |x| < \aj \t. 

Lax reflection, x < 0. For < 0, aj < 0, and > 0, we compute from Lemma [7] (followed by direct 
integration) 



"(*-*)>" _ _ 2 

\- ^TT-S H M ^-1/2,-1/2, ' 





(t - s) _1 e ~M(t- s ) (i _|_ s)- 1 /2 s -i/2 e -— dyds 



<Ct- x / 2 / (t-s)- 1 / 2 (l + s)- 1 / 2 e »* ds. 



In the case |x| > |a ■ |t, we write 

x — a ■ (t — s) — a, -^s = (x — a ■ t) + (a ■ — a;^r)s, 

for which there is no cancellation and we obtain an estimate of the form of a diffusion wave. We observe that 
for \aj /a^l > 1 the estimate arising from this integral will in general be broader than the diffusion wave we 
began with (we will get a constant greater than L). We can proceed, however, by dividing the integration 
over s into cases, s £ [0, t/C], C sufficiently large, for which we have control over the breadth of the diffusion 
wave, and s £ [t/C,t] for which we have exponential decay in t, which can be subsumed by our estimates. 
(For x > C\t, C\ sufficiently large, we clearly have exponential decay in both |x| and t; hence, we can take 
decay in t to give similar decay in |x|). 

For the case |x| < \a~\t, we first consider the case s £ [0,t/C], for which we write 

- a 7 - - - a 7 

x - a (t - s) - ai^s = (x - a t) + (a , - a l - J r)s, 

a, a. 
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and proceed as in the Lax convection case with s E [0,t/2]. For s E [t/C, t — t/C], we clearly have (t + \x\) 1 
decay, which can be subsumed, while for s E [t — t/C, t] we write 

x — a - (t — s) — a, -^s = (x — a, -^s) — a ■ (t — s), 
a k a k 

and proceed as in the Lax convection case with s E [t/2,i\. 

Lax transmission, x < 0. For a ; ~ < 0, > 0, and af > 0, we consider integrals of the form 

* f o (-• =t»-4*) a (B _ a _ s)2 

/ (t - s)-^ k — e-^il + s)^ 1 /2 s -i/2 e di—dyds 

J-oo 

a + 

(«-»+(t-.)-«r -Jr.) 2 



M((t-s) + (-^) 2 



In this case, for |x| > et, and fixed e > 0, we have exponential decay in both |x| and t, which can be subsumed 
into our estimates. In the case |a;| < et, for the case s E [0,t/2], we write 



x — af(t — s) — a, -^s = (x — aft) + (af — a, -^)s, 
a k a k 



and proceed as in the Lax convection case with s E [0, i/2], while for s E [t/2, t] we write 



x — af(t — s) — a^^s = (x — aV '-zrt) + (a^ - ^ — a t)(* — s )' 
a k a k a k 



and proceed as in the Lax convection case with s E [t/2, t] . 

Undercompressive convection, x < 0. For aj < 0, we consider integrals of the form 



(t - s)- l ' 2 e W^) e-" lal (l + s)- 1/2 s~ 1/2 e ^^dyds, 



t - s y 1/2 e~ 

la J-oo 

for which we observe the inequality 



e -v\y\ e - (v < Ce'^e-" 123 . 

We can now employ LemmaEJand proceed as in the Lax convection case, taking advantage of the integrability 
of e - '' 1 ^ in y and the integrability of e~' ,2S in s. Estimates in the undercompressive reflection case and 
undercompressive transmission case follow similarly. This concludes the analysis of the nonlinearity 9 2 . 
Nonlinearity ip 2 . We estimate convolutions of the form 

t M 

/ \G y (x, t - s; y)\{l + s^s^tfdyds. 

J-oo 

Lax convection, x < 0. For aj < 0, we consider integrals of the form 

(t _ a )-l e - l ""M(i-T n ' (1 + S )V2 s -l/2 



J a 1 s 

x (1 + \y\ + s)- 1 ^ + \y- ajsD^dyds. 

Clearly, the primary new element here is that we no longer have the precision afforded by Lemmas El and 
We proceed instead through balance estimates similar again to those that arise in the proof of Lemma El 
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We consider three subcases: (f ) a- < < a k , (2) a- < a k < 0, and (3) a k < a- < 0, beginning with the 
first. We begin by considering the interval |x| > |aj~|i, on which we write 

x — y — a k (4 — s) = (x — a^t) — (y — a^t) — a k (4 — s). (87) 

For s € [0, 4] and y G [a^s, 0], we have y — a^t > and consequently each term in 187(1 is negative, and we 
have no cancellation. Integrating in straightfoward fashion, then, we can conclude an estimate by C6{x,t). 
For |x| < |di"|£, we begin by writing 

x — y — a^(t — s) = (x — a^(t — s) — ajs) - (y - ajs), 

from which we see that either x — y — a k (t — s) is near x — a k (t — s) — aj s or 

\y-ajs\ > e\x-a^(t-s) -ajs\, 
for some e > 0. More precisely, we have the estimate 



(l + \y- a, s\)~ 1/2 e 



< C 



(»-»-«.- (t-s))" 

(»-»-a- (t-s)) 2 



(1 + \x - a k (t - s) - ajs\y 1/2 e 



+ (l + \y — a J s\) 1/2 e • < 

For the second piece of this estimate, we have 

I (t-s)~ 1 e Af(t- S ) J e ~ " ««-«) 

J a~ s 

x (1 + sf^s- 1 ' 2 ^ + | v | + s)- x (l + |y - ajslj-^yds. 

On the interval s S [0, 4/2], we integrate (1 + |y — ajs|) _1 in y, while on the interval s S [4/2, 4] we integrate 
the remaining Gaussian kernel. We obtain an estimate, then, by 

ft/2 (:c-a~ (t-s)-aT s) 2 

C / t- 1 (l + s) 1/2 S - 1/2 log(l + C's)e M '(*-=> ' ds 



n 



( x -"fc (t— ») — O ■ s) 2 

+ C (l + t)- 1 (t-s)- 1 / 2 e «'<*-> 3 cfe, 

Jt/2 

both of which can now be analyzed by the methods of the Lax convection case with nonlinearity 6 2 . 
The critical new estimate in this case is 

ft rQ (a; _ H _ a - (t „ s)) 2 

/ / {t-s^e {l + \x~a-(t- s)-ajs\)- 1 / 2 

Jo J a 1 s 

x (1 + s ) 1 / 2 s - 1 /2(i + \ y \ + s )-i(i + \ y _ arsD-Va^. 
For |x| > |a~|4 (though |x| < |o7/|t; we recall that the case \x\ > \aj\t has already been considered), we write 

x - a k (* - s ) - °7 S = (?- a j l ) - ( a k - a 7)(* ~ s )> 

for which we have no cancellation and straightfoward integration provides an estimate by (1 + 4)~ 1 / 2 (1 + 
Ix-a'tl)- 1 / 2 . For \z 
s e [0, t/2], we write 



r — a k t\) 1 / 2 . For \x\ < |o- |4, we divide the analysis into the cases s € [0,4/2] and s G [4/2,4]. In the case 



x - a fc (4 - s) - a 3 s = (x - a fc 4) - (a fe - ^ )s, 
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for which (observing that for x < and a k < 0, \x — a k t\ = \x\ + \a k \t), 

(l + \x- a-(t - s) - aj8\)-V*(l + \y\ + s)^ 2 
< C[(l + |^| +t)- 1 / 2 (l + |y| + s)- 1 / 2 

+ (1 + \x - a- k {t s) - aJ S \)-V\l + \y\ + \x\ + t)" 1/2 ' 
For the second estimate in lj89|) . we integrate (1 + \y — a~s|) -1 / 2 to find 

/ / (t-s)- x e (l + s^s^il + lyl + s)- 1 / 2 

x (1 + \y - ajs\)- 1/2 (l + \x-a^(t-s)- ajs\y 1/2 {l + \x\ + t)- 1/2 dyds 

r t/2 

< Ct- X (l + \x\ + t)- 1 ' 2 / (1 + sf^ds 
Jo 

<C<- 1/2 (l + |a;|+i)- 1/2 . 
The first estimate in <|89ll can be analyzed similarly. In the case s £ [t/2,t], we have 

x - a k (t - s) - ajs = (x - ajt) - (a k - aj)(t - s), 



for which we have 



(t- S )-^ 2 (l + \x-a k (t- S )-a- S \) 

<C\\X- a^r 1/2 (l + \x-a7(t- s)a;s\)- 1 / 2 



{t-s)- 1 ' 2 {l + \x~ajt\)- 1 ' 2 



For the first of these, we estimate 



-1/2 



(89) 



(t - s)- l ' 2 \x - ajt\-^ 2 (l + \ x -a^{t- s)ajs\ 

It/2 Jajs 

(zs-y-a~ (t-s)) 2 

x e MirrT) (i + s ) 1 /2 s -i/2( 1 + | y | + s y^(iyds 
^Cil + ty^x-ajt]-^ 2 

x [ (1 + \x - a k {t - s)aT s |)- 1 /2( 1 + s )i/2 s -i/2 ds 

Jt/2 

<C(l + t)-V 2 \x-aJt\-V 2 , 

and similarly for the second. We remark that the apparent blow-up at x = a~t is an artifact of the approach 
and can be removed by the observation that for \x — a*Jt\ < C\/t, we can proceed by alternative estimates 
to get decay of form 9(x, t). 

The remaining cases aj < a k < and a k < aj < follow similarly as in the Lax convection case for 2 , 
with the arguments augmented by balance estimates along the lines of JHHl- 

Lax reflection. For aj < 0, a k > 0, and a~ < 0, we consider the convolutions 



t r0 



(^--i-a-a. (t-s)r 



(t-s) _1 e ~M(t-s) 



-/a, s 



x (1 + s) 1 / 2 s- 1 / 2 (l + \y\ + s)"^! + \y - aTs^dyds. 



For \x\ > \a ■ \t, we write 



x 3 -y - a ■ (t - s) = (x - a , t) - (-^y - a ■ s). 
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Here, x — a - t < and -^y — a - s > 0, and hence we do not have cancellation, and we get an estimate by 
C9(x,t). For \x\ < \aj\t, we write 

a~ _ _ _a~ a~ 

x - -^y - a- (t - s) = (x - a j (t - s) - a l -±s) - -^(y - a ; s), 



from which we have the inequality 



{ X --2-y- a . (t-s)) 2 

(1 + \y - a- S \)- 1/2 e ^Mits) 



< C 



(1 + \x - aj(t - s) - a," — s|r 1/2 e 



Ik 

M(t-s) 



{x-a (t-«)-o, -±s) 2 (x- -J-y-aT (t-s)f 



+ {l + \y-ais\)-^e 



M'(t-s) 



M(t-s) 



We proceed now as in the Lax convection case for writing for s G [0, t/2] 



and for s E [t/2,t] 



x — a,j (t — s) — a ; -^s = (x — cij t) + (a,j — a l -^)s 



_a- a- a- 

x — cij (t — s) — a l -^s = (x — a l -^t) + (a ; — )(t — s). 



Lax transmission. For a ; < 0, a k > 0, and > 0, we consider the convolutions 



f f {t-s)-\ 

J J a, s 



(»— i7«-«J(t-»))' 

M(t-s) 



-rj\x\ 



x (1 + a) 1 /2 a -V2( 1 + | y | + s )-i(i + \ y _ a-sD^dyds. 

In the case |x| > et, some fixed e > 0, we have exponential decay in both \x\ and t, which can be subsumed. 
In the case |a;| < et, we proceed almost exactly as in the Lax reflection case for tpf. 
Undercompressive convection. For aj < 0, we consider the convolutions 



. y-q fc (t-s)) 

M(t-a) e -v\v\ 



I [ (t-s)-^e-- 

JO J a\ s 

x (1 + s) 1 /2 a -V2( 1 + \ y \ + s )-i(i + \ y _ a^sD^dyds. 
We observe here the inequality 

e -ri\v\(i + \ y _ ajs])- 1 < Cle-^e- 7128 + e'^l + s)" 1 

Integrating e~^ y \ (or e~^\y\), we now proceed similarly as in the analysis of the Lax convection case for 
Similarly, the undercompressive reflection and transmission estimates follow as in the Lax reflection and 
transmission esimates. This concludes the analysis for the nonlinearity tpf. 
Nonlinearity ip 2 - We estimate convolutions of the form 

f-t f-aZ s 



[ f * S \G v (x, t - a; y)\(l + s) 1 / 2 *" 1 / 2 ^?/, sfdyds. 

JO J-oo 
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Lax convection. For a 1 < 0, we consider convolutions of the form 



/ / (t-s)-h 

JO J -co 



(s-g-g-tt-s))" 
M(t-s) 



We write 
for which 



x (1 + s)V2 a -V2(i + \ y _ a - s \ + s l ' 2 )- s dyds. 
x-y-a k (t-s) = (x - %(t-s) - a^s) - (y-a^s), 



(x — y — a. (t-s)) 2 

(1 + \y - a^s\ + s 1/2 )" 3/2 e 



< C 



(1 + \x -a k (t-s)- a7s| + s 1/2 r 3/2 e 



(x-a^ (I-.)-! °) 2 .(«-„-»- (t-s)) 2 

+ (1 + \y - aj" s| + s 1/2 ) _3/2 e 5?7 T^> e 



M(t-s) 



For |x| > |a x |i, we write 

x - a k (t - s) - ajs = (x - a^t) - (a k - of)(t - s), 
for which there is no cancellation and integration gives an estimate by 



Ci(l + i)- 1 log(l + t)e- l '~£ t ty + C 2 (l + t)-V4(i + \ x _ a - t \ + t l ' 2 Y 



-3/2 



For \x\ < \a 1 \t, we consider only the case a k < 0. The case a k > is similar. We first consider the additional 
subcase \x\ < \a k \t, for which we write 

x - a k (t - s) - a^s = (x- a k t) + (a k - a7>. 

Observing that there is no cancellation between these terms, we can proceed exactly as in the case \x\ > jaj" \t 
to get an estimate by C(9(x, t) + tp 2 ). For the case \a k \t < \x\ < \aj\t, we divide the analysis into subcases 
s e [0,t/2] and s e [t/2,t]. For s e [0,i/2], we write 



a; - a fc (i - s) - a x s = {x - a k t) + (a k - a x )s, 



for which we have 



(1 + \x - a k (t - s) - a 1 s\ + s 1/2 ) 



-3/2 



< C 



(l + lx-a-tl+s 1 / 2 )- 



3/2 



(90) 



+ (1 + \x - a k (t - s) - aj>| +\x- a^ 1 / 2 + s 1 /2)-3/2 
For the first of the estimates, integrating (1 + \y — a{s\ + s 1//2 )~ 3 / 2 , we estimate 

ft/2 ra>i s 



/ / (t-s)- 1 (l + \x-a k (t-s)-a^s\+s 1/2 )- 3/2 e 

JO J-OO 

X (1 + S ) 1 /2 S -1/2 (1 + \ y _ a - s | + s V2)-3/2 dyds 
,t/2 

< Ct- 1 / (1 + \x - a k t\ + ,s 1 / 2 )- 3 / 2 (l + S ) 1 / 2 ,s- 1 / 2 (l + s 1 ' 2 )- 1 ' 2 ^ 



<c{i + t)- z i\i + \x- ak -t\y 1 / 2 . 

For s G [t/2,t], we write 

a; - a k (t — s) — a±s = (x - a^t) - (a k - a^)(t - s), 
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for which we have 

-3/2 



(t - s)-^(l + \x- a- k {t s) a- lS \ + s^ 2 ) 



< C 



\x~a^t\- 1/2 (l + \x-a-(t-s) -oJfs| + s 1/2 ) 



+ (t-s)-^ 2 {l + \x-a-t\ + S 1 / 2 ) 
For the first, integrating over the Gaussian kernel, we estimate 



3/2 

a 1 s\ + s x/ " I 

—3/2-1 



t fa, s 



(t - s)~ x (l + | at - a^(t - s) - a^s\ + s 1/2 )~ 3/2 e 

t/2 

X (1 + S ) 1 / 2 S - 1 /2(1 + \y _ ars | + S l/2)-3/2 d2/ds 

<C(l+t 1 /3)-3/2| a ._ art |-i/a 

' (1 + |z - a^(< - s) - a-s\ + ^/2)-3/2 (1 + s) l/2 s -i/2 ds 

t/2 

<C*(l + <r 1/2 (l + |:E-a^i|r 1/2 . 

The second estimate in ()90|l can be analyzed similarly. 

The Lax reflection and Lax transmission estimates follow similarly. Finally, for the undercompressive 
estimates in the case of nonlinearity ^f, we observe that for y G (— oo, a^[s] exponential \y\ decay yields 
exponential s decay, and the estimates follow in straightforward fashion. This concludes the analysis for the 
nonlinearity ip 2 and consequently for the nonlinearity 

Nonlinearity (1 + s) (9 + ipx + t/^). We consider convolutions of the form 

rt rO 

\G y {x, t-s;y)\(l + s)" 1 ^ + + fajdyds. 



Lax convection. For a,j < 0, we consider integrals of the form 

(»-t-»~ (t-s)) 2 ( H -a~s) 2 

(f - s) _1 e (1 + sy 1 s~ 1/2 e n dyds. 

We observe that this integral is better than the one analyzed in the Lax convection case with nonlinearity 
6 , except that the constant L in the diffusion kernel must now be kept. (In general, our balance estimates 
increase the size of this constant, see especially Lemma 0) According to Lemma we can estimate this 
integral by 

ct- 1 ' 2 [\t - s )- 1/2 (i + a )- l e - ^<«- 1 .)+^. cx-»>(«-)-»r') 3 ( fa. 



For | a: | > |a ■ \t and the case a- < a k , we write 

x — aj~(t— s) - ajs = (x — ajt) - {a k - aj)(t - s), 
for which there is no cancellation and integration yields an estimate by 

{l+t)- 1 / 2 e J ^M>t ) , 

where M < M' < L. For |sg| < \aj\t, we divide the analysis into the subcases s S [0, t/Ci] and s £ [t/Ci,t] 
for some constant C\ sufficiently large. In the event that s € [t/Ci,t], we integrate (t — s)^ 1 / 2 to obtain the 
estimate 

C(l+t)-\ 



28 



P. HOWARD and K. ZUMBRUN, Stability of Undcrcompressive Shock Profiles 



which can be subsumed. On the other hand if s E [0,i/Ci], wc observe that by choice of C\ the divisor 
M(t — s) + Ls can be kept as close to Mt as we like. Since M < L we can use this observation to recover the 
expected estimate. Otherwise, the analysis proceeds exactly as in the Lax convection case for nonlinearity 
9 2 . 

The Lax reflection and Lax transmission estimates for nonlinearity (1 + t)~ 1 6 follow from the Lax con- 
vection argument and Lemma The undercompressive estimates follow similarly. 
In the case of the nonlinearity (1 + s) _1 -0i, we observe that for y 6 [a^s, 0], 

(l + S )- 1 <C(l + |y| + ,s)- 1 , 

and so 

(l + s)-Vi<C^. 

Hence the convolution estimates for his nonlinearity follow from those for the nonlinearity ip 2 . Finally, 
estimates for the nonlinearity (l+s) -1 ^ are straightfoward. This concludes the analysis of the first estimates 
of Lemma [3J 

Excited term estimates. The remaining estimates of Lemma regard the excited terms e(y,t). For the 
integral 

/ / \e y t(y,t- s)\y(y,s)dyds, 

JO J -oo 

we have according to Remark 

\e yt {y,t)\<C{t- l +1 t- l l 2 e-^) £ 



(«+° fc t) 
e Mt 



a k >0 

We observe that these estimates correspond precisely with the Lax and undercompressive convection kernels 
with x = 0. In this way, we immediately obtain an estimate by 

C(0 + Vi+?fe)(M) <C x {i + t)- x , 
for some constant C\ . For the integral 

/■oo />0 

/ / K(y,+oo)|*(2/, s)dyds, 

JO J -oo 

we have from Remark [B] 

\e y (y,+oo)\ < C 7 e-"W 

where as usual 7 is 1 for undercompressive profiles and otherwise. For the nonlinearity (1 + s) 1 / 2 s _1 / 2 6' 2 , 
we observe the inequality 



e -v\v\ e 



(B-o. «)* 



< Ce-^e-^, 



for some fixed rji > and 772 > 0. We estimate 

roo rO 



e -m\v\ e -V2'(l + s)- 1 l 2 s- 1 / 2 dyds < C, 

by the integrability of e _?,1 ' a ' and e~ ms . Similarly, for the nonlinearity (1 + s) 1 / 2 s _1 / 2 -0 2 , we use the 
inequality 

e -v\v\(i + \ y _ ajsl)- 1 < c\e- 7>lM e- r>2S + ^"^(I + s)" 1 ] , 

from which the required estimate follows from the integrability of e~ ?n ' a ' and of (I + s)~ 2 . For the nonlinearity 
(I + s) 1 / 2 s _1 / 2 -0|, we have y 6 (— 00, aj~s], for which exponential \y\ decay yields exponential s decay and 
the estimate is immediate. For the nonlinearity (1 + s)~ 1 (9 + + 1/12), we observe the inequality 



e -nl»l/ 



' + +fo)(y,s) < Ce-^il + sy 1 . 
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We have, then an estimate by 



C / / e~ r ' lM (l + s)- 2 dyds 



For the integral 
we have from Remark 



J -oo 

/>oo 

< Ci / {1 + s)- 2 < C 2 . 
Jo 

t M 

\ \e y (y,t- s) - e y (y, +oo)\^(y, s)dyds, 

J-oo 



\e y (y,t)-e y (y,+oo)\<Ct- 1 ^ £ e 

a7>0 



Observing that our estimate on \e y (y, t) — e y (y, +oo)| takes the form of the Lax convection kernel multiplied 
by (t — s) 1 / 2 , we proceed exactly as there to determine an estimate by C(l + ty 1 / 2 . 
For the integral 

/>oo 

\e y (y,t - s)\V(y, s)dyds, 



we have from Remark 

\e y {y,t)\<Cr^ e " 



a, >0 



+ C 7 e-"lf I (errfn(^^) - errfn(^=^i)) . 



For the first, we can proceed simiarly as in the Lax convection estimates for x — 0. For the second, we 
observe the inequality 

e -"l»l*(y,») < Ce- mM {l + s)-^ 2 s-^ 2 . 

We estimate, then 

oo r o 

m\v\h j- ^-3/2 c -i/2 



e-^(l + s)-^ 2 s-^ 2 dyds 

) 

< c(i + ty 1 ' 2 / (i + syh-^ds < c'(i + ty 1 / 2 . 



This completes the proof of Lemma [3] 
I 

Proof of Lemma |4j Observing that the first three estimates of Lemma 0] are significantly less precise 
than the second three, and moreover can be established by similar methods, we begin with the fourth, 

ft rO 

/ \G y (x,t- s;y)\$ 2 (y, s)dyds, 

J—ac 

where 

$2(V,«) < Ce"^(l + s)- 3 / 2 . 
Convection estimate. We consider integrals of the form 

f [ (t - s )-W e - lX ~ V "<t-->' n2 e-"l"l(l + s)-^ 2 dyds. 

JO J-00 

In order to make use of the localization due to the term e - ''' 2 '', we use the inequality 



n?=7) e -v\y\ < c 

(»->-ar(t-.))" 



(*-<»-(»-.))•' 



g M'(t-s) e ~V\v\ 

M(t-s) e —n\y \ e -V2\x-a~ (t-s)\ 



(91) 
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For the first of these two estimates, we have the integral 

f f (t - s)- 1 /^" <X M'tt-s? e -" lal (1 + s)- 3 / 2 dyds. 

JO J -oo 

In the case |x| > |aj"|t, we have 

x-aZ{t-s) = (x-aTt) - (al - a7)t + als 
< (x-a7t) < 0. 

We have, then, an estimate by 

Cxt-We—T&r- / {l + s)- z ' 2 ds 
Jo 

Jt/2 

< C(l +t)- 1/2 e- <X . 

For |x| < \aT\t, we write 

x — aZ(t — s) = (x — aZt) + aZ s. 
In the event that aZ > 0, we have \x — aZt\ > c(\x\ + 1), so that 

(i + s y 3/2 e ~ ( M>( t -s ) )) < C (\x\ + ty 3/2 . 

In this case, we have the estimate 



C{\x\+t)- 3 ' 2 / (t-sy^ds^CQxl+t)- 1 . 



For a k < 0, we first consider the case \x\ > \a k \t, for which we again write 

x ~ a k (* ^ s ) = ( x - a k *) + a k s - 

Observing that there is no cancellation between these terms, we can conclude an estimate by C9(x, t) as in 
the case \x\ > \aT\t above. In the event \x\ < \a~Z\t, we integrate (for e > sufficiently small) 

f f (t - s)- 1/2 e~ <X e ^ lyl (1 + s)- 3 / 2 dyds 

JO J-oo 

pe\x—aZt\ ( x -a~t) 2 

^c^- 1 ' 2 e — ^(l + S y 3 / 2 ds 



Jo 

+ c 2 r ] / 2 / (i + s )- 3 ' 2 ds 

J e\x — a k t\ 

ft/ 2 (x-a~t) 2 

+ Czt- 1 ' 2 e ^(l + s )-3/2 rfs 

JC\x-a k t|At/2 



+ c 4 (i + t)- 3 / 2 f 

hi; 



(t a)-W 

t/2 



< Ci(l + ty 1 l 2 e J -^ L - + C a (l + <)- 1/2 (l + \x - a-t])- 1 ' 2 . 



For integration over the second estimate in lj91|l . we proceed similarly. 

In this case the analyses of the reflection estimate and the transmission estimate do not differ significantly 
from the analysis of the convection estimate. 
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Overcompressive excited estimates. For the integral 



t r 



\e t (y, t - s)|e~" lal (1 + s)- 3 / 2 dyds, 



J -oo 



we have an estimate by 



t r o 



Cl I (t - s)-V 2 e- (B+ M( t -.* )) e ~v\v\ (i + s )- 3 / 2 dyds, 

JO J -oo 



for which we observe the inequality 

(y + a.7 (t-s)) 2 



M(t-s) e -WI < C 



e -ml»l e -»?a(*-») 



We have, then 



t rO 

/ (t - .s)- 1/2 e- ??1 l y le~ T ' 2 ( f - s )(l + s)- z ' 2 dyds 

J-oo 

<C ie -T f [ t/2 (l + s)- 3 / 2 ds 
Jo 

+ C 2 (l + ty 3 / 2 f (t-sy^e-^-^ds 

Jt/2 

< c{i + ty 3/2 . 

\e(y,t- s) - e(y, +oo)|e" ? ' l! ' l (l + s)- 3/2 dyds, 



For the final integral 



we have an estimate by 

C f f errm( l2/| " fz_ s) ) e -^(l + s^dyds, 

Jo J-oo My/t - S 

for some constants M > 0, a > 0. In this case, we observe the inequality 

errfn( |y| ~ Q ji^ )e-^l < Ce^ 1 '^- 1 ' 2 ^, 
Myt — s 

from which the estimate follows as above. 

The remaining cases of Lemma 0] can be analyzed similarly. a 
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